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Abstract 

Some mathematical models for the estimation of the effects of CrylAb 
and CryfF Bt-maize exposure in the biodiversity are discussed. Novel 
results about these models are obtained and described in the note. The 
exact formula for the proportion of the population that suffers mortal- 
ity exposed to CrylAb pollen, underlining its dependence on the margin 
from the Bt crop edge, is derived. In addition, regarding CrylF pollen 
effects, it is proposed a procedure, using a probabilistic and statistical ap- 
proach, that computes the width of the non Bt-stripes used as mitigation 
measures. Finally, it has been derived a lower bound, using probabilistic 
consideration, on the species sensitivity of Lepidoptera. 

1 Introduction 

In the last years, the cultivation of genetically modified plants (GMP) is very 
widespread in the world, in particular in America and Asia. On the other side, 
the debate in the scientific community and in the public opinion about the 
GMP effects is becoming harder and harder [7j HJ [12]. It is our belief, that 
one of the causes of the debate depends on the lack of available mathematical 
models that allow to assess quantitatively the effects of the GMP cultivation 
on the biodiversity. To this purpose, recently appeared some approaches that 
discuss the effects of Bt-maize |H El E3 E] , and two mathematical models for 
estimating the effects of CrylAb and CrylF Bt-maize [8] [10], on non-target 
lepidoptera. To the best of our knowledge [9] , these latter models are the 
first and the unique mathematical models that assess quantitatively the effects 
of Bt-crop on the biodiversity. For this reason, we have analyzed in detail the 
models [5J [TO], obtaining novel results. In particular, regarding the CrylAb 
maize [8] has been derived the exact formula, and not the approximated as in 
[8], for the proportion of the population that suffers mortality P. The exact 



formula allow to enlight the effective dependence of P on the margin from the 
Bt crop edge D. 

Regarding CrylF model [ID] , we have derived a procedure, using a probabilistic 
and statistical approach, that allows to fix automatically the width of the non 
Bt-stripes used as mitigation measures. Moreover, we have shown that the 
values of width, proposed in 10J, does not implement the worst-case method, 
recommended by the directive of the European Community. Finally, it has 
been derived a lower bound, using probabilistic consideration, on the species 
sensitivity of the hypothethical Lepidoptera considered in [ID] . 
The work is organized as follows: In Section [2] the model for the estimation of 
CrylAb Bt-maize effects is discussed and the exact formula for the proportion 
of the population that suffers mortality P is derived; In Section [3] the model 
for the estimation of CrylF is analyzed and the procedure that computes the 
width of the non Bt-stripes is presented; finally, some conclusions are drawn in 
Section [4] 



2 Estimation of CrylAb Bt-maize effects 

Perry et al. [5] developed a mathematical model for estimating the effects of 
exposure of three non-target Lepidoptera (e.g., Inachis io, Vanessa atalanta, 
Plutella xylo stella) to Bt-maize pollen expressing the protein CrylAb. They 
computed the estimated proportion of the population of a non-target Lepi- 
doptera that suffers mortality, hereinafter denoted by P for convenience. P is 
given by 

yzvxa(2heh\J C + fDfi) 
~ (25eVC + fD) ' 

where the parameters h (within-crop mortality), x (physical effects), a (temporal 
coincidence), z (maize cropping), v (utilization rate), y (host plant in arable), 
e (host plant within-crop), / (host plant in margin), C (size of maize fields), 
D (width of margin) are experimentally measured or taken by the literature; 
whereas the parameter [i, depending on D, has to be properly computed. 
Following Perry et al., before computing /i it is necessary to recall the so-called 
margin mortality g(E) where E denotes the distance from the edge of the crop, 
that is given, for larvae of Inachis io and Vanessa atalanta, by 

whereas for larvae of Plutella xylostella is given by 
9 (E)= CM-0M9E) 

yy ' 55.33 + exp(-0.349i;) v ; 

That being said, Perry et al. declared that /x is obtained by averaging the value 
of g(E) over the margin. In practice, \x is obtained by the numerical integration 
of equations ([2| or ([3]), between the values E = and E = D. In this way, they 



2 



obtained an approximate estimate of [i and, consequently of the proportion of 
population that suffers mortality P. Besides, not having an exact formula that 
states fi in function of D, they cannot establish the effective dependence of the 
proportion P on the margin D. 



2.1 The exact proportion suffering mortality formula 

Now, we show that, differently of what claimed by Perry et al., an exact value 
for fi and hence P can be derivable. The parameter fi is given by 

D J HV D J 5 + cxpi-^E) 

where 7 and 5 are 0.359 and 33.25 for the larvae of Inachis io and Vanessa 
atalanta, whereas are 0.349 e 55.33 for the larvae of Plutella xylostella. 
The integral of the equation Q can be solved exactly: 

1_ f D cxp(- 7 £) 
Wo 5 + eM-lE) 
1 f D r S 



D J 1 6 + exp{--fE) 1 



1 f» SexpfrE) 



D J 1 Sexp( 1 E) + l' 
= i[£-~ln(l + <5exp( 7 £0)]£ 

= ^[D--\n(l + Sexp(-fD)) + -\n(l + S)} 

Plugging last formula in the equation ([T]) , we obtain the exact formula for com- 
puting the proportion of population that suffers mortality P 

p = yzvxa{25ehVC + f[D i In 1+S ffi D) ]} 
25eVC + fD 

The previous equation allows us to predict quantitatively, and not qualitatively 
as performed in Perry et al.'s work, the proportion of population that suffers 
mortality P at large distance D from the crop edge. We pass to study the 
dependence on D of the proportion P. When D = 0, i.e., when the non-target 
Lcpidoptera is on the crop edge, P becomes: 

P(0) = yzvxah. (7) 

When D goes to the infinity the equation ^ reduces to: 

yzvxa(25ehy/C + L m 1±£) x 



P(D)f*~ = J - —^O - . (8) 

V ' 25eVC + fD \Dj 
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This means that to the infinity P goes to as j~j that implies that even for large 
D, P is not negligible. 

3 Estimation of CrylF Bt-maize effects 

Perry et al. jTUj developed a mathematical model for estimating the effects of 
exposure of five hypothetical non-target Lepidoptera species to Bt-maize pollen 
expressing the protein CrylF. In the work the estimated proportion of the pop- 
ulation of a non-target Lepidoptera that suffers mortality, P, is given by the 
equation (1), that for our convenience, we recall: 

. . yzvxa(25ehy / C + fD/i) 

The previous equation combines: (a) small scale parameters, namely, e (host 
plant within-crop), / (host plant in margin), C (size of maize field), D (width of 
margin); (b) large scale parameters, namely, y (host plant in arable), z (maize 
cropping) , x (physical effects) , a (temporal coincidence) ; (c) mortality parame- 
ters, namely, h (within crop mortality) and /i (average mortality within a margin 
of any particular width D). Small and large scale parameters are taken from 
literature or experimentally measured. The mortality parameters, instead, have 
to be properly computed, hence in the following the focus will be on the way 
this is accomplished. 

Equation ([!]) is derived for a range of five hypothetical non-target Lepidoptera 
species rather than a specific one. To this aim, Perry et al. introduce a further 
mortality parameter, m, representing a range of species sensitivities for the hy- 
pothetical non-target Lepidoptera. Sensitivity is expressed by the LC50 values 
for maize 1507, i.e., the lethal concentration value that kills on average half of 
the larvae of the instar considered, measured in pollen grains per cm 2 . The 
parameter m affects the way h and fi are computed. Let us focus on a step-by- 
step mortality parameter derivation [10] . The starting point is a mortality-dose 
laboratory-derived bioassay relationship in which a logit-transformed probabil- 
ity of mortality, P, is regressed on a logarithmically transformed dose, d: 

logit(P) = a + 2.473 log 10 d. (9) 

Here is where the parameter m comes into play. The intercept a is determined 
by the sensitivity of the species to the CrylF protein, for which logit(P) = 0. 
Five sensitivity values, corresponding to five hypothetical species, are considered 
in 1 1 Oj and denoted as worst-case, extreme (m = 1.265), very high (m = 14.36), 
high (m — 163.2), above-average (m = 1853), and below-average (m = 21057). 
These mortality-dose relationships are then combined with a field-derived re- 
gression of logarithmically transformed dose, d, on distance E, from the nearest 
source of the pollen: 

log 10 d = 2.346 - 0.U5E, (10) 



4 



to derive a linear mortality-distance relationship for mortality of larvae in the 
margin, on the logit scale. So doing, for each species sensitivities the mortality- 
distance relationships are derived, from equations Q and (10), as 



logit{P) = A) - 0.3586£, (11) 

where (3q — a + 5.8017 and whose values, for each of the considered sensitivities, 
are shown in Table [TJ 



Sensitivity 


A) 


extreme 


5.5492 


very high 


2.9399 


high 


0.3297 


above-average 


-2.2798 


below-average 


-4.8901 



Table 1: Intercepts, A, for the linear mortality-distance relationships in equa- 



tion (111, corresponding to the considered species sensitivities. 



Taking the inverse of the logit function, we return to the natural scale thus 
obtaining the estimated probability of mortality g(E), for a larva at distance E 
into the margin from the nearest source of pollen at the edge of the field: 

yy ' y v ' 1 + cxp(/3 - 0.3586£) (3 + exp(-0.3586£:) v ' 



In equation (12), the corresponding values of f3 = exp(— A), for each species 
sensitivity, are listed in Table [2] The mortality parameters h and /j, are derived 
from g(E). 



Table 2: Values of j3 to compute g{E) in equation (12 1, corresponding to the 
considered species sensitivities. 



Sensitivity 




extreme 


0.0039 


very high 


0.0529 


high 


0.7191 


above-average 


9.7747 


below- average 


132.9669 
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Concretely, Perry et al. compute /i by numerically integrating g(E), in equation 
(12), between and D, as described in [5]. 

To estimate the probability of mortality, h, for the five hypothetical larvae within 
the -Bi-crop, it is necessary to consider that pollen deposition within a maize 
crop is 2.757 times that at the edge [5]. Therefore, 

h = 2.757 5 (0) = 2.757^. (13) 

The values obtained for h, for the considered range of sensitivities, are shown 
in Tabled 



Table 3: Within-crop mortality probabilities, h, corresponding to the considered 
species sensitivities. 



Sensitivity 


h 


extreme 


2.7463 


very high 


2.6185 


high 


1.6037 


above-average 


0.0928 


below-average 


0.0075 



To conclude the model overview, two further parameters have to be introduced, 
i.e., mitigation parameter and large scale exposure parameter. We describe only 
the former, since the latter does not affect our discussion. 

Perry et al. [10] considered mitigation measures, too. Indeed, the parameter w 
is the width in metres of the non- Bt maize strips which represents the simulated 
mitigation measures which are assumed to be planted around each of the four 
field edge. When there is mitigation (w > 0), similar calculations to those 
just described are used. However, mortality calculated for larvae in the margin 
must use an appropriate value of E, computed to allow for the fact that the 
Bi-maize is a distance w metres further away. The same reasoning is made for 
the mortality of larvae within the non-Bi-maize. 

3.1 Computation of the non-Bt maize strips width 

In this Section, we propose a procedure that fixes automatically the width of 
non Bt-maize strips. To this purpose, it is necessary to consider that the CrylF 
Bt-maize exposure model is based on the CrylAb Bt-maize one, previously de- 
scribed in Section [2j As a consequence of that, the model inherits from CrylAb 
Bt-maize exposure model, the proportion of population that suffers mortality 
P. Therefore, the analysis regarding the computation of the exact value of P, 
described in Section [2j is equally valid for the CrylF Bt-maize exposure model, 
too. This implies that P, given by the equation Q, asymptotically tends to 
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jj where D is the distance from the crop edge. That being said, we propose 
a method to assess automatically, starting from the expression of the equation 
(|8j), the width w of the non-Bt maize strips. We propose to use the 3a rule [TT] . 
borrowed by the statistical and probability theory, to derive a fixed value for 
the width w . Let be P(0), see equation Q, the proportion of population that 
suffers mortality on the crop edge, the width w is given by the value of D such 
that: 

where 77 is the probability that a zero-mean normal distribution with standard 
deviation a assume values in [— 3<r, 3a], i.e. ~ 0.9973. Since P(D) « O(j^), for 
large D, we fix: 

w = {D: I«(l-ry)} (15) 

Hence D can be fixed to 1/(1 — 0.9973), i.e. ~ 370m. If we use this D value 
to fix the width w of the non-Bt maize strips we have the consequence that 
the Bt-maize field, that we assume for simplicity squared, must have a size 
larger than twice the value of D, i.e., 740m. For instance, we consider the 
example described in [TO] with the unique difference that the size of the field 
is 774m (the double of the one considered by PerrjQ. It's easy to show that 
the area of the field that can be devoted to the Bt-maize cultivation is a square 
of size 34m (see Figure [l} , whereas the rest of the field cannot be used since it 
should be destined to mitigation field. The example shows that the adoption 
of mitigation measures implies strong constraints. Firstly, field size should be 
larger than 740m; secondly, only a limited field area can be devoted to the Si- 
maize cultivation. We conclude the analysis showing how the values proposed 
by Perry et al. can be viewed in our statistical and probabilistic framework. In 
particular if it considers the most severe value proposed by Perry et al., namely 



24m, after the rule (15) application it obtains for rj a value of ~ 0.0416 that 
corresponds to apply approximatively a 2a rule (2.03<r [3] for the precision) . 
This choice cannot be absolutely considered to fulfill the worst case scenario, 
recommended in the Directive 2001/18/EC of the European Community [2]. 

3.2 Derivation of a lower bound for species sensitivity 

Now we examine the within crop mortality parameter, h, and its assumed val- 
ues listed in Table [3] of the section [3j Three out of five h values, corresponding 
to the extreme, very high and high sensitivities, violate the definition of prob- 
ability. Recall that the probability alway takes on values in between zero and 
one |llj . while the above mentioned h values are greater than one. The joint 
assumptions that pollen deposition within a maize crop is 2.757 times that at 
the edge, and the considered range for the species sensitivities do not appear 
to fit the theoretical soundness of the proposed model. The right side of equa- 



tion (13) should be < 1, hence equation (13) defines a probability if and only 



x the original size of the field considered by Perry et al., 387m, cannot be considered since 
it is lower than 740m. 
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Figure 1: Crop Field. The parameters A, p, w stand for the size of squared 
crop field, the size of i?i-maize cultivated field, the width of non-Bt strips, 
respectively. 



if P > 1.757. Therefore, more restrictive assumptions on the sensitivity should 
be made and/or further argumentation should be given concerning with the as- 
sumption that pollen deposition within a maize crop is 2.757 times that at the 
edge. Obviously, the incorrect values of within-crop mortality h, for extreme, 
very high and high sensitivities are propagated through the model affecting the 
correctness of the proportion of population that suffers mortality P computed 
in equation ^ and, consequently, of the entire model. However, assuming the 
above mentioned Perry et al.'s issues are reasoneable, a lower bound for the 
values of sensitivity m can be added to the Perry et al's model in order to get 
a correct probability value for the parameter h. 

We pass to derive the lower bound. By applying the definition of probability, 
the following relations must hold: 

2 757 

Kl» - — - < 1 =>■ > 1.757, 
1 + p - 

and, since P = exp(— Pq), we get 

exp(-A)) > 1-757 po < - In 1.757 => p < -0.5636. 



It comes, from equation (11), that Po — a + 5.8017 and therefore, a < —6.3653. 
Now, from equation ^ we get: 

logit(P) - 2.473 log 10 d < -6.353 (16) 
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and setting logit(P) = 0, it follows: 



2.473 log 10 d 

logio d 
d 



> 



> 



< 



2.5739 



-6.3653 
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,2.5739 



374.89. 



(17) 



Therefore, since at LC50, m = d (see supplementary material SI and S3 in 
[TU]). the lower bound m > 374.89 must hold in order to h to be a probability. 

4 Conclusion 

In this work, some mathematical models for the estimation of the effects of 
CrylAb and CrylF Bt-maize exposure on non-target Lepidoptera have been 
discussed, deriving novel results. Firstly, it is obtained the exact formula for the 
proportion of the population that suffers mortality exposed to CrylAb pollen, 
studying its dependence on the margin from the Bt-crop edge. Besides, regard- 
ing CrylF pollen effects, a procedure was proposed to fix automatically the 
width of the non Si-stripes used as mitigation measures. It was also shown 
that the adoption of mitigation measures, that take into account the worst-case 
scenario recommended in the directive 2001/18/EC by European Community, 
implies strong constraints on the Bt-maize cultivation. Firstly, field size should 
be larger than 740m; secondly, only a limited field area can be devoted to the 
Bt-maize cultivation. Furthermore, on the basis of probabilistic considerations, 
it has been derived a lower bound, on the species sensitivity of Lepidoptera. 
We hope that the derived novel results about mathematical models for the es- 
timation of the effects of CrylAb and CrylF Bt-maize exposure on non-target 
Lepidoptera can contribute to the debate in the scientific community about the 
GMP effects. 
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